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Abstract 

Standard methods in non-linear analysis are used to indicate that 
there exists a parabolic branching of solutions of the Lichnerowicz- 
York equation with an unsealed source. We also apply these methods 
to the extended conformal thin sandwich formulation and, by assum- 
ing that the linearised system develops a kernel solution for sufficiently 
large initial data, we reproduce the parabolic solution curves for the 
conformal factor, lapse and shift found numerically by Pfeiffer and 
York. 



1 Introduction 



With the onset of gravitational wave astronomy approaching, it is of crucial 
importance to have constructed suitably realistic theoretical and numerical 
models of the space-time structure of gravitational wave sources. This de- 
mands that we study the initial value problem for Einstein's gravitational 
field equations. These form a complicated quasi-linear system of partial dif- 
ferential equations. They naturally split into six evolution equations and four 
constraint equations on the initial data. It is the constraint equations that 
this work will focus on. 

*E-mail : darraghmw@gmail.com 
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This work was initiated to try to explain the intriguing non-uniqueness 
results found numerically by Pfeiffer and York in [TJ. They found a parabolic 
curve of regular solutions for each of the five unknowns in the extended 
conformal thin sandwich formulation of the Einstein constraints. This sys- 
tem has enjoyed much support amongst numerical relativists but their results 
show the existence of two regular solutions for the entire range of wave ampli- 
tude considered. We present a simple local derivation of these non-uniqueness 
results. 

We first review the important features of the conformal method for solving 
the constraints. The four constraint equations of general relativity are 

R(g) - fC'Kij + K 2 = 16np (1) 

Vi{K ij - g ij K) = 87cJ j (2) 

where is the extrinsic curvature of the space-like initial data slice and 
K = tr-gK l i and p and J- 7 must satisfy conservation equations. In the various 
conformal formulations an initial 3-metric gij is chosen which is conformally 
related to the physical solution of the constraints g~ij 

9ij = <fr 4 9ij- (3) 

We decompose as K l i = A^ + \g l ^K and define the trace-free conformal 
extrinsic curvature tensor A % i according to 

A ij = <j)- 10 A ij (4) 

which leads to 

A ij Aij = (j>~ 12 A ij Aij, (5) 

where barred objects are defined with respect to the physical metric g and 
K = K. This property is used in the transformation of the Hamiltonian 
constraint ([1]). Together with the transformation for the scalar curvature 
given by 

R{g)<f = R{g)<P - 8V 2 (6) 
it yields the Lichnerowicz-York (LY) equation for 0: 

V 2 - ^ = —A^- 7 + ^ - 2^. (7) 
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Under the conformal transformations ©-(jl]) a symmetric trace-free tensor 
satisfies 

ViA" = <f>- W ViA ij , (8) 
which simplifies the momentum constraint (J2]), 

ViA ij - -<fS7 j K = 87r^ 10 J j . (9) 
3 

We see immediately that the vacuum (i.e. p = J 3 ' — 0) constraint equations 
decouple if we have constant mean curvature, K=const. 

The original conformal transverse traceless (CTT) formulation by York 
and his later conformal thin sandwich (CTS) formulation have identical ex- 
istence and uniqueness properties (see [2]). They differ in their construction 
of the tensor . The CTT method relies upon tensor splittings to con- 
struct a TT tensor. The CTS formulation is much simpler and bypasses 
these complications. There A l i takes the form 

A" = ^- U " (10) 

where is the trace-free part of the time derivative of the conformal metric, 
N is the conformal lapse related to the physical lapse N by N = N<ft 6 , (3 l 
has a natural interpretation as the shift vector for the spatial slice, and L is 
the conformal killing operator defined as LJf y ' = V*X J + V^X 1 — |<7 y 'VfcX fc . 
The initial data for the four equation system (j7|), flUJ) is (g^, U^, K, N). 

The extended conformal thin sandwich (XCTS) system [I], [3] extends 
the CTS system by allowing the slicing condition to be propagated (K is 
now initial data so that the evolution equation for K becomes a constraint 
equation). This is highly desirable since it is natural to choose K — at 
quasi-equilibrium (i.e. when C/ u = 0), for example, while it is unclear what 
choice for the conformal lapse in the standard CTS formulation should be 
chosen for this situation. However, this condition couples the lapse fixing 
equation K to the four constraint equations and yields a far more complex 
system with no constant mean curvature decoupling. The initial data is now 
in the Lagrangian form (g^, Uij; K,K). The five vacuum XCTS equations 
are 



v2 Rm^_«w£ +«l 

8 32x 2 YT 

32 X + ^2~ 



V\ - » = ^ + ^ -^(K- ^K) (12) 
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n 2Ar y 3 v ; 

where a(/3) = (L/3 y - U i j)(hf3 ii - U ij ) and x = N<p 7 . This extended system 
shares the conformal covariance properties of the original system by construc- 
tion since the new equation it couples to is the physical lapse fixing equation. 
In the analysis that follows we restrict our attention to asymptotically flat 
initial data. 

The non-uniqueness results in [lj should not be confused with the non- 
uniqueness that results from choosing trivial initial data in the standard 
CTT/CTS formulations. In this case we need only solve ([7]) which is simply 
V 2 = 0. We get a unique regular solution (flat space = 1) or any number 
of singular moment of time symmetry Schwarzschild solutions by linearity. 
However, Pfeiffer and York found two regular solutions for each non-zero 
wave amplitude considered and their results suggest that the limiting case 
of trivial initial data is identical to the CTT/CTS scenario (as it should be 
given that the two systems are essentially equivalent for trivial initial data). 
We only consider regular solutions in this work. 

In the next section we review the basic uniqueness properties of the max- 
imal (K=0) CTT/CTS formulation of the constraints. In section three we 
focus on the LY equation with unsealed sources. We assume that a (critical) 
solution exists whose linearisation has a kernel and show using Lyapunov- 
Schmidt theory that a parabolic branching of solutions. In section four we 
consider the non-uniqueness results for the XCTS system. We show that with 
the assumption of a ID kernel, the full XCTS system exhibits a parabolic 
branching in all five variables, exactly as was found numerically in pQ. We 
conclude with a discussion of the possible implications of these results for 
evolutions. 



2 The Linear System and Bifurcations 

It is important to understand why the branching features described in [1] 
and outlined below do not occur in the maximal CTT/CTS formulation. 
It is known that the LY equation admits unique solutions [4J away from 
trivial initial data. For the maximal (i.e. decoupled) CTT/CTS systems the 
existence and uniqueness of solutions of the LY equation is determined by 
a variational inequality for the initial data. A theorem by Cantor and Brill 
in [5J, later corrected by Maxwell in \6\, states that on an asymptotically 
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Euclidean manifold there exists a positive solution of the LY equation if and 
only if V / € the following inequality holds 

mff e cg° j|7ij| > (14) 



where R is the scalar curvature of the conformal metric and the volume 
element and inner product are with respect to this metric. No such theorem 
is known for the XCTS system. When studying the p0 5 equation and later 
the XCTS system we will emphasise that both systems arise from variation 
of a non-convex action so that solutions would correspond to saddle points 
and are not expected to be unique. 

In non-linear systems the existence of a bifurcating branch of solutions is 
indicated by the non-invertibility of the linearized system. Thus in order to 
explain the non-uniqueness of solutions in [1] we need to find a background 
solution whose linearisation is not invertible. 

Given a nonlinear scalar equation Au = f(x,u), such as fl7|) with K=0, 
our first task is to check if it is linearisation stable, see [8]. If f jU < then 
the linearisation has the 'wrong sign' for use of the maximum principle to 
show local uniqueness of solutions. To illustrate this we recall the following 
result proven in [9]: 

On a suitably smooth asymptotically flat background the operator 

(V 2 -/): H k>5 -> ff fc -2,*- a (15) 

where / > and f £ H 2j s with 5 Q < —2 (so that f is continuous and 
falls off faster than r~ 2 , these function spaces are defined in the appendix) 
is an isomorphism if — 1 < 5 < and 2 < k < 4. If we linearise the 
maximal CTT/CTS vacuum LY equation about a solution it reduces to this 
isomorphic form, so that the implicit function theorem gives us complete 
neighbourhoods of nearby solutions ('linearisation stability'). This provides 
a qualitative explanation for the uniqueness results found by Pfeiffer and 
York for the maximal 4-equation CTS system (see Fig. 1 in pQ). 

In the next section we consider the unsealed source model studied by 
York in [TU] which is not linearisation stable, it's linearisation has the 'wrong 
sign' for application of the maximum principle so that a kernel solution may 
exist for sufficiently large initial data (the function f in ffT5l) is negative). We 
demonstrate the existence of a parabolic solution curve similar to the ones 
found in [TJ in a small neighbourhood of such a critical solution. 
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3 Non-Uniqueness for Unsealed Sources 



In [10] York considered the case of moment of time symmetry conformally 
flat initial data with positive energy density p (a more general analysis of 
unsealed sources including existence proofs may be found in [UJ). Then the 
LY equation (JTj) reads 

V 2 + 27rp0 5 = (16) 

where <fi — > 1 at spatial infinity. He noted that this equation is not lineari- 
sation stable because the linearisation has the 'wrong sign' for use of the 
maximum principle to prove local uniqueness. 

We now study the form of this non-uniqueness using Lyapunov-Schmidt 
(LS) theory. When a linear operator B has a kernel we are no longer able 
to invert the equation BV=h(x), with h(x) given, unless fV*h = where 
V* is the kernel of the adjoint problem i.e. V* G KerB* = coker B. The 
LS theory may be regarded as a local extension of these ideas to non-linear 
problems. In this method we remove the kernel of the linearisation from 
the domain and project the source terms, now non-linear combinations of 
the unknown V, onto the image of the linearisation. In this way a bijective 
operator B is defined in (1221) . The implicit function theorem gives a unique 
solution to the mapping between the modified spaces from which we recon- 
struct a solution to the original problem from a Taylor expansion of a real 
valued function (the LS equation). 

We briefly review the LS theory here following closely [12] (see also [13] ) 
whilst keeping in mind the immediate application to eqn ffTol) . We work in 
weighted Sobolev spaces (defined in the appendix) to ensure that the integrals 
below are finite, because the Fredholm properties of mappings between these 
spaces are well known (see [9], [2D]) and because the implicit function theorem 
is easily defined for such spaces. We write the nonlinear equation in functional 
form F(X, A) = where F is a smooth function of the unknown X and 
a parameter A. F defines a map between Hilbert spaces (defined in the 
appendix) F : X x R — > Y. To apply the implicit function theorem F must 
be at least C l and this follows from the multiplication properties of functions 
in these spaces. For an example of the implicit function theorem applied to 
the CTS formalism with between these function spaces see [UJ. Once 
a solution to the restricted problem is obtained via the implicit function 
theorem the remainder of our analysis is formal. 

Our analysis is perturbative and we will assume in this section and the 
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next that there exists a critical solution X c corresponding to parameter value 
A = A c . By this we mean that when we linearise F about (X c , A c ) we find 
that the linear operator B := DxF(X c , A c ) has a kernel i.e. there is a non- 
zero function Vq such that BVq = 0. We require that B be fredholm of index 
zero (i.e. Dim Ker B=Dim coker B, see appendix) and focus our attention 
on the case when the kernel of B is ID. We perturb A c by A = A c — e and 
look for solutions to F(X, A)=0 of the form X = X c + V where V is a small 
perturbation. 

Clearly F = is equivalent to 

BV = R(V,e) (17) 

where R(V, e) = BV — F(X, A). For equation ([TBI) B is self-adjoint and (for 
functions with falloff satisfying 5 G (—1,0), see appendix) we have V = V$ 
where Vq spans the kernel of B*, the formal adjoint of B. Given Vq G H^ g we 
define z G C£°, for simplicity, so that j Vqz = 1. We then split the domain 
of B as 

XcH ktS 3V = £V + u. (18) 

If we define the parameter £ = j Vzdv then j uzdv = 0. Thus X is split 
parallel and perpendicular to z by this choice of parameter. We take k > 4 
so that if V G H^s then we also have V G C 2 (see the appendix). 
Similarly we split the range of B as 

^C# fc _ 2 , 5 _ 2 3R = dz + w (19) 

and taking d = J RVodv means J wVodv = and so the Range is split parallel 
and perpendicular to Vq. 

This splitting allows us to define the following projection operators 

P : = X -> X, PX = £V (20) 
Q : = Y -> Y, QY = dz. (21) 

We now apply separately the projections 1 — Q and Q to ([17]) to obtain 

Bu = R(£V + u, e) (22) 

= QR(£V + u,e) } (23) 

where the operator B : (1 — P)X — ► (1 — Q)Y is now bijective. This means 
that the linearisation of f[2"2l . Bu = 0, only has the zero solution. The implicit 
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function theorem then implies that there exists a unique small solution u = 
«(£, e) G C 2 to (I2"2"j) with u(0, 0) = and u^(0, 0) = since there are no terms 
in (122]) linear in £. 

We now substitute this result into (123]) which yields the LS equation (see 

QR(£V o + u(£V ,e),e)=0. (24) 

With u(£Vq, e) a known function given by the implicit function theorem, this 
should be viewed as a real valued equation. It gives gives a relation between 
the parameters £ = £(e). 

Since Q is a projection operator we may write the LS equation as 

d(£,e)=/ i2(£F + u(£Vo,e),eW« = 0. (25) 
i.e. d = in (fT9D . 

At this stage, from the splitting (fT8l) . we have a solution to ffTTj) of the 
form 

\Z = £(e)y + M(£,e) (26) 

with the particular form of £(e) determined by (1251) . We note that for each 
value of the perturbation e, £(e) traces a curve along which the splittings 
above are valid. If D\F((p c , A c )e ^ for some e then we know that the set of 
zeros of d(e), d _1 (0), form a smooth submanifold (see [TB]). 

We now proceed with a formal argument to determine £(e). Expanding 
the known C 2 function u as a Taylor series in the dependent parameter £ 
and the independent parameter e we have 

u(£, e) = eu e (0, 0) + 0(2) := eu* + 0(2) (27) 

since ug(0,0) ^ is incompatible with (|2"2"|) due to the absence of linear 
terms in V in R (where subscripts denote partial derivatives and 0(2) denotes 
second order terms in £, e). Now our solution to (|T7|) is 

V = Z(e)V + eu* + 0(2). (28) 

The LS equation may now be written in the form 

L 20 + £eLn + eL i + ... = (29) 

where L mn = J VoR mn dv and R mn denotes the mth order term in £ and 
the nth order term in e resulting from substitution of the solution (1281) into 
R(V,e). This fixes £ = £(e). 
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In [12] Vainberg and Trenogin prove that the small solutions (i.e. where 
£(e) — > as e — > 0) of the LS equation are in 1-1 correspondence with the 
small solutions of F(X, A) = 0. 

We now show that branching of the type found for the XCTS system in 
PQ is actually a generic property of solutions of the standard LY equation 
with an unsealed source (this was shown recently in the special case of the 
constant density star in [21]). We examine the local behavior of the solutions 
to this equation at a critical point of the linearisation as the unsealed source 
term is varied. We multiply p in ([16]) by a positive parameter A and seek a 
continuous family of solutions to the LY equation 



on a fixed flat background. If p is a suitably smooth compactly supported 
function and we seek solutions such that — 1 G H k § where 5 G (—1, 0) 
then the integral relations above will be finite. 

For A = the only regular solution satisfying the boundary conditions is 
= 1. As we increase A we find that the linearised homogenous equation 



has only the trivial solution Vo = 0. The implicit function theorem then tells 
us that fl30|) has a smooth sequence of solutions, 0(A), with 0(0) = 1 and 
the maximum principle tells us that 0(A) > 1. We now assume that there 
exists a critical solution C at A = A c , whose linearisation has a kernel Vo, 
corresponding to the lowest eigenstate of the Schrodinger-type equation (T5TT) 
. We know that this eigenstate is unique up to scaling and has no nodes, so 
we may take Vo > 0. 

We now expand (|30|) about this critical value where the linearisation, 
f !3ip . has a kernel. Taking 

:= C + V with C > and A := A c — e we find 



V 2 + 27rpA0 5 = 0, ^ 1 as r ^ oo 



(30) 



V V + lO7rpA0 4 V o = 



(31) 



V 2 (0 C + V) + 2np(\ c - e)(0 c + V) 5 = 



(32) 



which gives 



BV 



(V 2 + 1OttA c 0^)V (33) 
2irep<f>l - 2OvrA c p0^'V 2 + Wnp^eV + 7 other terms. (34) 
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The LS equation (1251) for this problem is 

/ Vq (2vr p<pl + lOnp^eV - 20tt \ c p<p 3 c V 2 + ...) dv = 0. (35) 

If we substitute the solution V in the form (1251) into (|35p we obtain the 
coefficients L mn : 

L i = / V (27rp<^)tfo > (36) 

L n = I V (107rp<f) 4 c Vo-407r\ c p<f) 3 c VoU*)dv (37) 

L 20 = [ V (-207r\ c p<j)lV 2 )dv < 0, (38) 

and so on. 

By choosing £ and e small enough we may truncate the LS branching 
equation (1291) at any order. In particular we can write 

£%o + £e|£n|+e£oi « 0. (39) 

Solving the quadratic equation for £ we find, 

i 

£ = ±L^ e y+ (Vi). (40) 

where — > as e — > 0. This tells us that to lowest order and provided 
Loi 7^ that we may ignore the contribution of Lu. 

Therefore, in a small neighbourhood of the critical solution <ft c we find 
that as the parameter e is varied that the conformal factor traces a parabola 
in the solution space 

= C ± (jl^ 2 V + O(V~e). (41) 

If L i = then 2irpcf) 5 c G Image(B) and a qualitatively different situation 
arises. The LS equation then tells us that there is more than one branch 
of solutions passing through A = A c provided not all 7^ 0. In this case 
the zeros of the bifurcation equation ( |35l) do not form a smooth submanifold. 
Note that this is not the behaviour observed in [1] where only a single smooth 
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(parabolic) curve of solutions is found. This phenomenon was observed re- 
cently in a numerical study of non-unique solutions to (IT51) corresponding to 
a constant density star, [13]. Further details can be found in [T3] . 

As noted above, specialising (Tl6]) to the constant density star, upper and 
lower branches of solutions were found in [21] and an explicit form for the 
critical and kernel solutions were given. Generalising from p = const and 
compactly supported to just p > and compactly supported, only the lower 
branch of solutions was given and a parabolic approach to the critical solution 
was found in [21]. The existence of an upper branch was conjectured and 
is easily determined using the LS methods outlined above which give the 
required ±-y/e behavior in a small neighbourhood of the critical point. 

In [10] York notes that by specifying a conformal energy density p re- 
lated to the physical energy density by p = p<f)~ s where s > 5, we transform 
the linearisation of fl30|) into an isomorphic form yielding unique solutions. 
Therefore non-uniqueness results from a poor choice of conformal scaling in 
this case. We argue in the next section that the coupling of the K evolution 
equation to the four constraint equations in the XCTS system similarly in- 
troduces an undesirable scaling of the variables leading to non-uniqueness as 
above. 

4 The XCTS System 

As in the case of unsealed sources above, to apply the LS method we must 
first assume the existence of a critical solution about which the linearisation 
of this system has a (vector) kernel. 

We note here that the XCTS system arises from the variation of an action 
given by 

A V g 12 AV g V J 

(42) 

We may also include scaled versions of the energy density and the current 
density. This allows us to examine the existence and uniqueness properties 
of the XCTS system from a variational perspective. 

First we investigate the necessary conditions for the functional (j42|) to 
have a minimum. We must have that the functional is bounded below on an 
appropriately defined domain of admissible functions II which will be some 
Sobolev space H k ^ subject to the constraints iV > 0, (p > 0. 
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Provided such a domain has been defined, we proceed to the calculation 
of the first and second variations of X. If we define J(s) = X{u + sv), 
where v G II, a stationary value of X is given by J'(0) = where ' denotes 
differentiation with respect to s. A simple calculation reveals that the XCTS 
system corresponds to a stationary point of this functional. 

It is straightforward to show that the second variation of (|42p is not of 
definite sign. The following inequality is used in [T5] as the basic assumption 
on the integrand F to apply direct methods in the calculus of variations; 

F p , p ,(x, M ,V«)£4>0 (43) 

for all rank one matrices ^ where i=l-5, a— 1-3 and P % a denote the usual 
derivatives in the Euler-Lagrange equations. This condition, called 'strong 
ellipticity' in [J5], is the requirement that F(x,u,p) be convex with respect 
to p (p = Vu G M 5x3 ). The XCTS Lagrangian (j42p fails this criterion; it is 
convex in L/3 but, due to the mixed term V%V0, F is not convex in V</> or 
in Vx- 

This lack of convexity means we cannot expect stationary points, if they 
exist, to be unique. An analogous situation occurs frequently in non-linear 
elasticity where non-convex functionals are necessary in order to model buck- 
ling equilibrium configurations of materials which are known physically to be 
non-unique, i.e. an input stress can lead to numerous buckled states. This 
argument and the analogy to the unsealed source equation ([IB]) lends support 
to our assumption below that the linearised XCTS system develops a kernel 
for sufficiently large initial data. 

4.1 Non- uniqueness in the XCTS system 

In this section we assume that the XCTS system (fTT|) - (fT3"|) has a critical solu- 
tion Ac_and that the linearisation about this solution has a one dimensional 
kernel V q. To apply the LS method we must also check that the formally 
adjoint system has a kernel of equal dimension. We know from [5] that for 
sufficiently smooth initial data with very general falloff conditions at spatial 
infinity that the kernel of this linearised system is finite dimensional and 
that it has a closed range. We show in the appendix that the linear system 
is actually Fredholm with an index of zero (so that dim Ker=Dim coKer) 
between suitably defined Sobolev spaces. Assuming now that the kernel is 
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one dimensional, we have satisfied the requirements to implement the LS 
theory outlined above. 

We introduce a parameter A to allow us to continuously vary the confor- 
mal initial data and in the conformal background(in [lj the parameter 
'A' was used which corresponded to the amplitude of a Teukolsky wave in the 
conformal background). The 1-parameter family of initial data considered in 
[ljis 

9ij = $ij + ^hij (44) 
Uij = Xhij. (45) 

where hij is the tracefree part of the metric perturbation (these tensors are 
given explicitly in [I]). This corresponds to a gravitational wave perturbation 
of flat space with a gaussian wave profile. The perturbation hy decays expo- 
nentially with distance so the metric is asymptotically flat. (Our construction 
is valid for systems with much weaker power-law falloff). Furthermore, due 
to the fast falloff, we know that this conformal metric will not contribute to 
the ADM energy of the physical solution. 

In ( !TT]) - (JI3"]) we denote the dependence of the initial data on A by 

V ^ V A , UV ^ Xlpi, R ^ R(X). 

When A = the XCTS equations (fTTl) - (lT3"l) decouple due to the non-existence 
of conformal killing vectors that vanish at infinity (where we have the con- 
straint on the conformal lapse that iV > 0), and we obtain flat space as the 
unique regular solution. Applying the implicit function theorem then gives 
a local curve of solutions parameterised by A. We assume that there exists 
a critical solution X c = (0 C , Xc, PI) occurring at A = X c . At this point we 
apply the LS method to continue the curve through A c . 

Below we will, as in the previous example, perturb the system about this 
critical solution. It is convenient to absorb the critical shift vector into U 1 ^ 
by a gauge transformation. We perturb the XCTS system (fTT]) - (jT3"]) at the 
critical solution X c = (0 C , Xc, 0, 0, 0) according to 

A = A c - e 

(j) = C + 01 

X = Xc + Xi 

p l = o + p[. 



13 



We expand the background scalar curvature in a Taylor series about A c 
so that 

R(X) = R(X C ) - eR'(\ c ) + ... (46) 
The linear terms in the expansion give the following inhomogenous system 



A^--R(X c)h + ^^-2U.m-(2^-7^)U.U) = e^ c + eT 
A X1 -^(A C ) X1 + 7^(2^ 1 -(6|-|)^) = e^W 

\ Xc Xc <Pc Xc J 

where terms arising from variation of the connection are given the generic 
symbol T. Due to the form of the conformal initial data we are considering 
we know that these connection terms decay exponentially with distance. 

The XCTS system (fTIl)-(II3I) will be referred to here as F(X, A) = 0. As 
before, F describes a mapping between Hilbert spaces F : X x R — > Y. We 
assume the existence of a ID kernel V o and define the dual object Z £ 
so that J R3 V o. Z dv = 1 where the inner product denotes multiplication of a 
five component row vector with a 5 component column vector and the volume 
element is taken with respect to the conformal metric g^. Then we may split 
the domain X as five copies of the one defined in the scalar case (|T8|) and 
similarly the range corresponds to five copies of (JT9l) . Likewise our relation 
for £ in the scalar case, namely £ = J Vzdv, becomes 

f = / V.~Zdv. (47) 

</R 3 

The critical solution satisfies F(X C} A c ) = and V o satisfies 

bv = d x f(x c , X c )V = 0. 

The LS machinery developed for the scalar equation ( |T6l) naturally generalises 
to systems of elliptic equations. As before we have that F = is equivalent 
to 

BV ■= D X F{X C ,X C )V = R(V,e). (48) 



14 



Following the same procedure as (I22p - (l23p . restricting the domain and 
range of B by applying the projections 1-Q and Q to R gives us 



B~u = R(£V + Tt,e) 
QR{£V + Tt,e) = 



(49) 
(50) 



where Q denotes the projection onto Z as before. The implicit function 
theorem is also valid for nonlinear systems (see for example [EE]). The lin- 
earisation of (l4"9l now reads B~u = 0. Since B is now an isomorphism, we 
obtain a unique small solution ~u(£,e) to (1411 with ^(0,0) = "??g(0,0) = 0. 
We substitute this into (!50|) to obtain the orthogonality relation that defines 
the curve £ = £(e). It is worth pointing out that, since the linear system is 
fredholm with index zero and the kernel is assumed to be one dimensional, we 
still have just one equation, (1501) . for one unknown £ (though the numerical 
value of £ now depends on all five unknowns). 

At this point we have a (vector) solution of the form 



just as before (see 



= £(e)V + e^* + O(2), 
. Equation (T3U|) reads 



(51) 
(52) 



v*.(R)dv = / (y;*, v*\ v*?\ V*? 2 , v* P3 ) 



( R* \ 

R x 
R^ 

\ R^ J 



dv (53) 



where R? corresponds to the first order in e background term, D\F(X C , A c )e, 
and nonlinear terms arising on the RHS of the Hamiltonian constraint. Note 
that in the appendix we have shown that 



/t/*<£ t/*x T/*/31 T/*/32 1ta*/33\ 



(V x V lp V &1 V 132 



Likewise the LS coefficients are given by 



La 



Vr>* .RiAv. 



(54) 
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where Rij denotes the ith order term in £ and the jth order term in e resulting 
from substitution of the solution f[5"Tj) into R . 

The source terms on the RHS of (HHD are of the form 



R) 



( 



\ 



-leR'(X c ) + eT + T 1 (<p, X ,L(3) + .. 
-leR'f(X c )+eT + T 2 ( ( j ) ,x,L(3) + . 

eT + T 3 ( ( f ) ,x,m + ... 

eT + T 4 (<f) } x,L(3) + ... 

eT + T 5 (cj ) ,x,L(3) + ... 



\ 



(55) 



/ R* \ 

Rx 
R^ 
R& 
\ Rfc ) 

where Tj represents the quadratic combinations of all variables arising 
from the perturbation. From the form of the solution (15TT) we see that the Tj 
terms yield terms quadratic in £ and e and mixed terms proportional to £ e. 

Clearly 

1 



-01 



K + V x )-R'(X c )-V .T)dv. 



(56) 



In the scalar model ffTB]) we worked on a fixed flat background so that the 
connection wasn't varied. We also knew that the first eigenfunction of the 
linearised equation had no nodes. When dealing with a coupled system of 
equations we lose this property. To obtain the results of [1] we must assume 
that L 01 7^ 0, which means that only a single curve of solutions passes through 
the critical point (<i _1 (0) is a smooth submanifold) . 

The next possibly non-zero terms are L\\ and L2o- If Lqi ^ and L20 7^ 
then we may truncate our series at quadratic order and ignore the contribu- 
tion of L\\. The complexity of the combinations of the quadratic terms in 
( 155]) do not yield easily to analysis. We prefer to emphasise that it is ex- 
tremely unlikely that our choice of initial data could lead to the cancellation 
of all terms at this order. Furthermore, since the kernel solution is removed 
by the LS method we know that the operator B is an isomorphism which 
implies that all the Ly's are finite for small £ and e. We conclude that 

£% + eL 01 « 0. (57) 
In principle the numbers L01, L20 can be determined for the choice of 
initial data that yields the kernel solution V (c.f. (1361) ). However, by 
choosing the sign of e so that < we know that there is a parabolic 
branching of all five variables 



X 



X c ± 



\L ie\ 
\L 



20 1 



V, 



(58) 
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(c.f. equation (7) in pQ). The parabola in pQ corresponds to the case where 
< so that e must be positive (i.e. there are no solutions with A > A c ). 
To recap, we list the assumptions that were made to derive fl58|) . The 
first assumption was that for sufficiently large A the curve of exact solutions 
(found numerically in [1]) reaches a critical solution (an exact solution whose 
linearisation has a 1-D kernel V o) and that none of the components of V o 
vanishes identically. If a component of was identically zero then the 
lowest order term in the expansion fl5~Tl) of this component would not vary 
like ±v^ i n the vicinity of A c , contrary to the results in pQ. 

Secondly, we assumed that L 01 ^ and L 2 o ^ 0. The very general form 
of the integral (1561) and the corresponding integral for L 2 o suggests that this 
is by far the likeliest outcome. This guarantees a parabolic curve such as 
( 158]) . The solution curves found numerically by Pfeiffer and York correspond 
to ^ < 0. 



5 Implications for evolutions 

We now consider the likely implications of these results for evolutions of the 
constraints. Recall that the Bianchi identities tell us that if the constraints 
are satisfied on the initial slice then they are satisfied on all slices of the 
foliation. However in numerical evolution schemes constraint violation off the 
initial data hypersurface remains a serious problem. To control constraint 
errors one may choose to solve the constraints on each slice of the foliation 
(constrained evolution). 

The Bianchi identities also make some of the Einstein equations redun- 
dant. Given a solution of the constraints (g = 4> 4 g,K), we may choose to 
evolve the conformal factor according to the evolution equation for the met- 
ric. For example, evolving initial data from the XCTS formulation naturally 
gives 

(d t - (5%)\og4> = h-N^K + V#) (59) 
o 

where N is the conformal lapse and derivatives are with respect to the con- 
formal metric. Whereas (|59|) selects just one solution of the constraints to 
evolve, the constrained evolution scheme could possibly lead to a solution 
jumping between branches during an evolution. 

From equation (l59i) we see that if we choose maximal slicing with zero 
shift then = and so a constrained evolution of this data should yield the 
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initial conformal factor throughout the evolution. This provides a possible 
test of whether a constrained evolution has caused to stray from the correct 
branch. 

The example f|T6|) serves to illustrate the need to choose an appropriate 
scaling of the extrinsic curvature. The standard scaling of the LY formulation 
in 01]) not only simplifies the momentum constraint but also removes the 
linearisation instability. Many axisymmetric evolution schemes (see [T7], |l8j ) 
do not scale the momentum. This need not cause problems on the initial slice 
if moment of time symmetry data is chosen. However, constrained evolution 
of this initial data could prove problematic as later time slices are susceptible 
to linearisation instability and non-uniqueness as noted in [19J. Also, the 
conformal scaling in the Hamiltonian constraint in the BSSN formulation is 
analogous to ffl6|) . This is not necessarily a problem if an initial data set is 
transformed into this formalism and then evolved using a free evolution such 
as (PJ. 

6 Conclusion 

In [1] the solution curves for initial data in CTS and XCTS were compared. 
As expected the CTS system had unique solutions for A < X cr u and no 
solutions for A > \ cr it- With identical initial data (p^, Uij), but with K = 
replacing N=l the solution curve for the XCTS system is parabolic and turns 
back upon itself at the (much smaller) critical wave amplitude. We have 
shown that this is consistent with the system developing a kernel solution of 
its linearisation and being dominated by a quadratic non-linearity there. 

It is worth emphasising that all the non-uniqueness results in this work 
are local in nature, i.e. the multiplicity of solutions is confined to a small 
neighbourhood of the critical solution. We have argued that the parabolic 
branching behaviour in [1] is the simplest example of branching phenomena 
that can arise in non-linear elliptic PDEs. However, the global nature of the 
parabolic branches in p] is surprising. 

The non-uniqueness properties described above are generic in the sense 
that they do not depend on the form of the (Fredholm) linearised operator 
B nor on the form of the forcing terms. If B has a ID kernel and there are 
quadratic nonlinearities then a parabolic solution curve is to be expected 
(provided L i, L 2 o 0). With higher order non-linearities or with multiple 
parameters Aj more complicated solution curves arise. It is therefore natural 
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to expect non-uniqueness (in a neighbourhood of the critical solution) in any 
non-linear elliptic system with a non-trivial kernel. 

Whilst our results using LS theory mirror the numerical results found in 
[1], it is important to remember that this non- uniqueness is removed when 
the wave amplitude A picks up a spatial dependence when we map from the 
conformal to the physical space. In the original paper it was noted that in 
the physical space the wave amplitude was not A but rather 4 A and that 
using instead Amax0 4 as our wave amplitude parameter leads to a solution 
curve with no turn-around point, solutions are seen to be unique for this 
"parameter" (see fig. 4 in [1]). We may draw a parallel to the scalar equation 
(I16p for which York noted that choosing an appropriate conformal scaling for 
p leads to linearisation stability and unique solutions but at the considerable 
cost of losing control over a physical quantity, the energy density p. 

The XCTS system has been very successful in modelling Black holes us- 
ing the puncture technique and trivial initial data corresponding to quasi- 
equilibrium i.e. (U^ = K = K = 0) for which the system considerably 
simplifies (see, for example, [22] where the value of the shift at the "punc- 
ture" was specified to prevent it vanishing everywhere). 

When we studied the XCTS system in this work we considered only 
Dirichlet boundary conditions on an asymptotically flat manifold containing 
gravitational waves. When modelling black holes an alternative to the punc- 
ture technique is to excise a region corresponding to an apparent horizon. 
Much analytical work has been done to guarantee unique solutions satisfying 
apparent horizon boundary conditions in the standard CTT formalism. A 
recent numerical study of the generalisation of these methods to the XCTS 
system revealed many difficulties with well-posedness (see [23] and references 
therein) as should be expected from the results outlined here. These results 
add weight to our reservations regarding the possible dangers of constrained 
evolution using ill-posed formulations of the constraints. 

When modelling data with Uij ^ 0, or when using apparent horizon 
boundary conditions it is clear that the XCTS system becomes much more 
complicated and it may prove more beneficial to revert to the standard and 
much simpler four equation CTS formalism. 
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Appendix: Function spaces and the linearised 
system 

We now define the function spaces used in the text. Our analysis has been 
restricted to asymptotically flat manifolds (M 3 ,^-) where g — 5 satisfies the 
falloff conditions below. 

We use weighted Sobolev spaces of tensors with norm defined as (see [20J) 



where a = (l+r 2 )a , r = \x\ is the Euclidean distance function and the volume 
form and covariant derivatives are with respect to the Euclidean metric 5. 

The conformal metric considered in [1] was a perturbation to flat space 
consisting of a quadrupole gravitational wave with a gaussian profile. The 
metric therefore had an exponential falloff with distance. When we confor- 
mally map to a solution of the constraints we know that this background 
metric will make no contribution to the ADM energy. 

Through the conformal method we aim to construct a Riemannian asymp- 
totically flat 3 metric g where 



and we seek solutions — 1, 1 — x, (3 l G H k ^ where k > 4 and 5 G (—1, 0). 
k is the number of times a tensor is weakly differentiable. Note that with 
k > 4 the metric is C 2 (we lose 3/2 degrees of differentiability in passing from 
weak to strong differentiability). If v G H k> s with k > 3/2 then \v (x) \ = o{r- 5 ) 
as r tends to infinity, see [2D] . 

The linearised XCTS system is a mapping 




(60) 



5ij G H k ,s 
Uij G Hk-i t s-i 



(61) 
(62) 



B : H, 




(63) 
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and was given after equation (l4"B"i) above. 

The formally adjoint system is defined for all u G H^^ according to 



<Bu,v>= J vBu = J uB*v (64) 

where inner products denote a product of five component row and column 
vectors, the volume form is with respect to the background metric and v G 
Hk-2 -1-5- Thus B* is a mapping 



B* : Iffc_2-i-<5 — > ,-3-s- (65) 



given by 



A0i - ii?(A c )0! + 7-^2- f 16E/.L& + (-6^ + —)[/.[/ ] = (66) 



32xo \ Xo' 
32 Xo \ ' X " ~Xo' 



Axi - ^(A c )xi + 77^0 [-16C/.LA + (7^-2^)U.u) = (67) 



L ^ + ^'(^-7^)l = 0. (68) 
Xo 8xo l Xo 0o ^ ' 

We now outline the Fredholm properties of B. Linear elliptic systems on 
asymptotically Euclidean manifolds of the form Bu = *Y^k=o CLkD k u were 
studied in [9]. Clearly the linearised XCTS system is elliptic. Our initial 
data satisfies hypothesis one of that work regarding smoothness and falloff. 
They proved that such a system has a finite dimensional kernel and a closed 
range. This says that B is semi-Fredholm so that the domain of B splits 
as s = kerB + W where W is closed and B is injective on W. However, 
this is not good enough for the application we have in mind. We need B to 
be Fredholm with an index of zero (so that the dimension of the cokernel is 
equal to that of the kernel of B). 

The domain of B is contained in H^^ such that 5 G (—1, 0). This implies 
that — 1 — 5 G (—1, 0) so that the domain of B* is the same as that of B. 
From a theorem in [9] we know that each system has a finite dimensional 
kernel and a closed range. At first sight it is not clear that the kernels of 
B and B* are related. However, inspection of the systems reveals that if 
V = (0i,Xl>/3i) G ker B then V* Q := (xij^i, §/?i) G ker B* and so 
the systems are identical under this relabelling of variables. Therefore the 
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systems have kernels of equal finite dimension so that B is Fredholm with an 
index of zero. This is the basic requirement for the use of the LS methods 
outlined above. If the linearised system is Fredholm of index zero then we 
know that the LS relation QR=0 gives n equations in n unknowns £j where n 
is the dimension of the kernel. If Dim Ker is 1 then we have a single equation 
for £ in terms of e. 
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